Mon. Not. R. Astron. Soc. 000, [|-?? (2001) Printed 1 February 2008 (MN 1*1^ style file vl.4) 



How old is the Universe ? 

Setting new constraints on the age of the Universe 

Ignacio Ferreras, Alessandro Melchiorri and Joseph Silk 

Nuclear & Astrophysics Lab. 1 Keble Road, Oxford OXl 3RH, United Kingdom 



o 
o 

(N 



(N 
(N 



> 

oo 
m 

o 



o 



Draft version 1 February 2008 



ABSTRACT 

There are three independent techniques for determining the age of the universe: via 
cosmochronology of long-Uved radioactive nuclei, via stellar modelling and population 
synthesis of the oldest stellar populations, and, most recently, via the precision cos- 
mology that has become feasible with the mapping of the acoustic peaks in the cosmic 
microwave background. We demonstrate that all three methods give completely con- 
sistent results, and enable us to set rigorous bounds on the maximum and minimum 
ages that are allowed for the universe. We present new constraints on the age of the 
universe by performing a multiband colour analysis of bright cluster ellipticals over a 
large redshift range (0.3 < z < 0.9), which allows us to infer the ages of their stellar 
populations over a wide range of possible formation redshifts and metallicities. Apply- 
ing a conservative prior to Hubble's constant of Hq = 70 ± 15 km s~^ Mpc^^, we find 
the age of the universe to be 13.2^2 Gyr (Ic), in agreement both with the estimates 
from type la supernovae, as well as with the latest uranium decay estimates, which 
yield an age for the Milky Way of 12.5 ± 3 Gyr. If we combine the results from cluster 
ellipticals with the analysis of the angular power spectrum of the cosmic microwave 
background and with the observations of type la supernovae at high redshift, we find a 
similar age: 13.4^j^'Q Gyr. Without the assumption of any priors, universes older than 
18 Gyr are ruled out by the data at the 90% confidence level. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: evolution — cosmology: 
Cosmic Microwave Background, anisotropy, power spectrum 



1 INTRODUCTION 

Discrepancies between age determinations, such as for globu- 
lar clusters and from the Hubble constant, have long plagued 
cosmology. The situation has changed dramatically in the 
past three years, however, as the Hubble constant propo- 
nents (Mould et al. 2000) have converged on a value which 
we conservatively take to be of Ho = 70±15 km s~^ Mpc~^, 
and because of several other developments. Type la super- 
nova measurements have provided strong evidence for a cos- 
mological constant (Perlmutter et al. 1999; Riess et al. 1998), 
and mapping of the acoustic peaks in the cosmic microwave 
background has fixed the cosmological parameters with un- 
precedented precision (Netterfield et al. 2001). Finally, de- 
tection of uranium in an old Population II star has pro- 
vided a direct nuclear chronometer for the age of our galaxy 
(Cayrel et al. 2001). Any one of these measurements may 
be suspect, but the remarkable concordance that we find 
enables to show here by a combined likelihood anaysis that 
combination of these constraints provides, for the first time, 
rigorous upper and lower bounds on the age of the universe. 

Consider first stellar model determinations of the age of 



the universe. The small scatter found in the absolute lumi- 
nosity of the brightest cluster galaxies (BCGs) over a wide 
range of redshifts motivated their use as standard candles 
to determine cosmological parameters (Gunn & Oke 1975). 
However, this analysis was based on the assumption that 
this type of galaxy should not undergo significant evolution 
in luminosity with lookback time. Tinsley (1976) showed 
that the brightening of main sequence stars poses a major 
hurdle in the use of BCGs as standard candles. Predicting 
the luminosity evolution is a rather challenging endeavour 
since it strongly depends both on the star formation his- 
tory as well as on the dynamical history. An analysis of 
bright ellipticals in a large sample of clusters observed in 
the near-infrared (Aragon-Salamanca, Baugh & Kauffmann 
1998) concluded that the stellar mass in BCGs over a large 
redshift range (0 < 2 < 1) has evolved by a factor between 
2 and 4 depending on the cosmology, in agreement with 
the predictions of hierarchical clustering scenarios of struc- 
ture formation (Kauffmann & Chariot 1998). On the other 
hand, an analysis of the colours of bright cluster ellipticals 
is only dependent on their star formation history. The small 
scatter found in the colour-magnitude relation of cluster el- 
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lipticals (Bower, Lucey & Ellis 1992; Stanford, Eisenhardt 
& Dickinson 1998) hints at old stellar populations, formed 
at redshifts zf ^ 3. The use of galaxy colours as a "cos- 
mic clock" is nevertheless a challenging task especially due 
to the age-metallicity degeneracy, which causes age effects 
to be mimicked by a range of metal abundances (Worthey 
1994). 

We focus on the brightest cluster ellipticals — which 
are assumed to have a simple star formation history as ex- 
plained below — and use stellar population synthesis models 
(Bruzual & Chariot, in preparation, hereafter B&C) in or- 
der to infer the age of the stars in these galaxies. We have 
compared this technique both with the analysis of the latest 
measurements of the angular power spectrum of the cosmic 
microwave background radiation (CMBR) observed by the 
BOOMERanG collaboration (de Bernardis et al. 2000, 2001) 
as well as with the observations of high redshift type la su- 
pernovae (Perlmutter et al. 1999; Riess et al. 1998) in order 
to estimate the age of the universe. We also compare these 
ages with the ages derived for globular clusters (Salaris & 
Weiss 1998) and the oldest halo stars (Cayrel et al. 2001), 
to which one has to add an age for the Milky Way that must 
correspond to the time elapsed between the Big Bang and 
formation at a redshift of at least 2, and more conservatively 
5 or even 10. 

We find that four completely independent age deter- 
minations — namely CMBR, galaxy colours for clusters at 
0^1, stellar evolution applied to old globular clusters, and 
radioactive isotope dating of old stars — lead to a consistent 
result. These age probes represent different combinations of 
Big Bang parameters {Hq, Sim, ^a), galaxy evolution pa- 
rameters (star formation rate history, initial stellar mass 
function), stellar evolution parameters (stellar mass, com- 
position and mixing length) , and nucleochronology (half-life 
of '^^^U), respectively. In this paper we add the spectropho- 
tometric study of bright cluster ellipticals to the growing 
list of cosmological probes used to determine the age of the 
Universe (Lineweaver 1999; Primack 2000). 



2 USING BRIGHT CLUSTER ELLIPTICALS 
AS COSMIC CLOCKS 

We use the sample of Stanford, Eisenhardt & Dickinson 
(1998), which comprises 17 clusters over a large redshift 
range (0.3 < z < 0.9). The sample was extracted on the ba- 
sis of available imaging with the Wide Field and Planetary 
Camera 2 on board the Hubble Space Telescope — for mor- 
phological classification purposes — from a larger sample of 
46 clusters drawn from a variety of optical. X-ray and radio- 
selected clusters. Each cluster was imaged in near-infrared J, 
H and K bands as well as two optical passbands, which were 
chosen as a function of redshift to straddle the 4000A break 
in the galaxy rest-frame (i.e. roughly mapping rest-frame U 
and V bands). For each cluster we select the three bright- 
est early-type systems which fall on the colour-magnitude 
relation and choose the reddest one, always taking care not 
to select an outlier. For each galaxy we compared the ob- 
served colours with the predictions of a grid of simple stellar 
populations with different ages and metallicities, from the 
latest models of B&C. The comparison was performed using 
a tsst applied between the four observed colours for each 



of the 17 clusters (c„,i; n = {1 • ■ ■ 17}, i = {1 • ■ ■ 4}), and 
the predictions from the population synthesis models (Cn,i), 
namely: 

Xssp = Z^lli Ylt=i (^"■^ ~ '^"•0 /"'^+ n\ 
+ (Model()7 - l^),=o - 1.6) Vc^o- 

The colour scatter ai is chosen to be ±0.2 magnitudes for all 
four colours, and comprises the effect of photometric error 
bars, uncertainties in the modelling of stellar populations as 
well as colour scatter in bright ellipticals. In order to further 
tighten the allowed region of parameter space, we have ap- 
plied a constraint at zero redshift, using high precision pho- 
tometry of the Coma and Virgo clusters (Bower et al. 1992), 
adding the second term shown in equation (1). Hence, for a 
given age-metallicity point, we evolve the system to zero red- 
shift (with a given cosmology) and then compare its U — V 
colour with the observed {U — V)z=o ~ 1-6 with a scatter 
ao = 0.1 mag. 

Figure 1 shows the contours at 68, 90 and 95% confi- 
dence levels in an age-metallicity diagram for the brightest 
ellipticals in cluster C10024-I-16 (z — 0.39). The degeneracy 
between age and metallicity is readily shown, so that old 
populations with low metal abundances give similar results 
to young stars with a higher metallicity. In the figure, we 
compare the result of B&C with the population synthesis 
models of Yi et al. (1999) who take special care in adding 
the contribution from core helium burning stars — horizon- 
tal branch (HB) stars — and their progeny. HB stars seem to 
be the most plausible candidate to explain the UV upturn 
(A - 1500A) in elliptical galaxies (O'ConneU 1999). How- 
ever, the agreement between both models shown in Figure 1 
is expected if we consider that most of the contribution of 
the light in the spectral range considered (i.e. between U 
and K bands) from elliptical galaxies comes from the main 
sequence and red giant branch populations, which are better 
understood than HB stars (Bruzual & Chariot 1993; Yi et 
al. 1999). 

We can impose a further constraint on the metallicity 
for bright cluster ellipticals. In these massive systems, the 
deep gravitational wells prevent gas ejected from supernovae 
from being thrown out of the galaxy. This feedback mecha- 
nism is presumed to be more effective in low mass spheroids 
(Ferreras & Silk 2000a) and is the basis for the correla- 
tion between mass and metallicity (Larson 1974) . A detailed 
analysis of metal abundances in several Fornax cluster early- 
type galaxies (Kuntschner 2000) shows that [Fe/H] ^ 0.0 for 
ellipticals with high velocity dispersion (log cr 2.2). The 
standard simple closed-box model is a good approximation 
to the chemical enrichment in bright ellipticals (Pagel 1997). 
In this model the evolution of the average stellar metallicity 
is: 



V 1 - /iJ 

where /i is the gas mass fraction contributing to star for- 
mation and p is the stellar yield, i.e. the mass fraction of 
elements other than helium generated in stars and weighted 
by the initial mass function (IMF) . In bright ellipticals there 
is no significant ongoing star formation, which means /i 
and so Z p <^ O.SZq (where Zq = 0.02 is the solar metal- 
licity). Hence, we can further constrain the upper bound to 
the age of the stellar populations by imposing a lower limit 



© 2001 RAS, MNRAS 000, |-?? 



How old is the Universe ? 3 



to the metallicity at [Fe/H] ^ — 0.1, roughly corresponding 
to Z ~ 0.8^0, shown as a horizontal line in Figure 1. 

Finally, the complete sample of 17 clusters can be com- 
bined by assuming a given cosmology {Ho,Qrn,^A) which 
enables us to translate redshifts into ages. A grid of models 
was run as follows: The cosmology was explored by choos- 
ing 0.085 < n-rn < 1.085 and < < 1., both in steps 
of 0.085; and 0.25 < ho < 0.95 in steps of 0.05 (with 
ho defined as Ho /km s~^Mpc~^). The stellar populations 
were parametrized by the metallicity, chosen in the range 
0.8 < Z/Zq < 1.9 in steps of 0.1 and the formation red- 
shift ZF = {2,2.25,2.50,2.75,3,3.25,3.50,3.75,4,6,8,10}. 
We ran two such models for two different initial mass func- 
tions: Salpeter (1955) and Scalo (1986), using B&C and a 
third one using the models of Yi et al. (1999) for a Salpeter 
IMF. 

The metallicity is not the only factor that could modify 
age estimates from broadband photometry. The reddening 
caused by dust may also cause overestimates of stellar ages. 
However, dust is not expected to play a significant role in 
cluster ellipticals, whose gaseous component is too hot to 
allow significant amounts of dust to be present over large 
timescales. The small scatter found in the colour-magnitude 
relation of cluster ellipticals would require a tight conspir- 
acy between age, metallicity and dust, in order to keep the 
scatter as low as observed. Furthermore, rest-frame near- 
ultraviolet photometry of Abell 851 {z = 0.41) shows that 
the brightest ellipticals are not redder than the dustless 
predictions of simple stellar populations (Ferreras & Silk 
2000b). 



3 AGE FROM COSMOLOGY: METHOD 

In the standard inflationary framework, the cosmic mi- 
crowave background power spectrum depends essentially on 
3 cosmological parameters (Efstathiou & Bond 1999): the 
physical matter density in baryons [cOb = f2i,/i^), the overall 
physical matter density {uJm ~ ^mh^) and the parameter 
R = ^uj^juJh^iy) where /(y) is sinh(j/), j/, sin(j/) for open, 
flat and closed models respectively and where 



y 



,1/2 



da 



(3) 



with oife = /i^(l — Vim — J^a). Cosmological models with the 
same values of the parameters ojm , '^'t and R will have nearly- 
identical power spectra on degree and subdegree angular 
scales. Furthermore, the age of the universe is given by 



to = 9.8Gyr 



ada 



(4) 



Under the assumption of a flat universe (i.e. ujk = Q 
and R ~ const) as the recent CMBR measurements seem 
to suggest, it is easy to show that nearly degenerate models 
with the same oja, L^m and lo^ will have also similar ages. 
Thus, in principle, a measurement of the CMBR spectrum 
can be extremely helpful in the determination of the age of 
the universe in flat cosmologies. The curvature itself, how- 
ever, due to the dependence of 7? from Oa and fim, cannot 
be constrained by CMBR measurements alone better than 
10 - 20%. 



This introduces a limitation in the use of the CMBR 
spectrum for producing independent strong constraints on 
the age. However, the integrated Sachs- Wolfe effect on large 
angular scales and the assumption of mild external pri- 
ors on the various cosmological parameters can break the 
above degeneracies and reduce the error in the age estima- 
tion from CMBR. In what follows we will put constraints 
on the age of the universe by comparing the recent CMBR 
data obtained from the BOOMERanG experiment with a 
database of models with cosmological parameters sampled 
as described in the previous section. We also vary the spec- 
tral index of the primordial density perturbations within the 
range = 0.50, 1.50, the optical depth Tc = 0.0, ...,0.3, 
and we rescale the fluctuation amplitude by a pre-factor 
taken as a free parameter. 

The theoretical models are computed using the publicly 
available CMBFAST program (Seljak & Zaldarriaga 1996) and 
are compared with the BOOMERanG-98 and COBE resuhs. 
We include the COBE data using Lloyd Knox's RADPack 
packages. The power spectra from these experiments are es- 
timated in 19 and 24 bins respectively, spanning the range 
2 < I < 1050. For the BOOMERanG-98 the spectrum we 
assign a flat shape, l{i + 1)Ci/2ti = Cb- 

Following de Bernardis et al. we approximate the signal 
Cb inside the bin as a Gaussian variable. The likelihood for a 
given cosmological model is then defined by L = c'^cmbrJ'^ 
with 



2 _ V^/^th 

XCMBR — / j\^B 



c^bYIoI, 



(5) 



where C^g (C^^) is the theoretical (experimental) band 
power, and ob is the quoted error bar. We consider a 10% 
calibration error for the BOOMERanG-98 experiment by 
adding a gaussian term xLi — (1-0 ~ ^cai)^/(0-24)^ and by 
finding the value of A^ai that for a given cosmological mod- 
els maximizes the likelihood. We also marginalize over the 
beam uncertainty (1.4') and we found that the removal of 
the last 3 bins — which are more likely to be affected by 
systematics — does not change the results of our analysis. 
We multiply L by our chosen priors and attribute a likeli- 
hood to each age in the 1 — 30 Gyr range by finding the 
'nuisance' parameters that maximise it. We then define our 
central values and Icr for the age from the 16%, 50% and 
84% integrals of L over age. 



4 DISCUSSION 

The results are summarized in Table 1 using external priors 
based on theoretical restrictions as well as on recent astro- 
nomical observations. The finite volume of parameter space 
sampled imposes further implicit constraints on the age of 
the universe. However, the large range of parameters ex- 
plored implies this implicit constraint is rather weak, as can 
be seen from the last entry in Table 1 (constant likelihood 
across the parameter space without any priors), for which 
the estimated age of the universe at a 68% confidence level 
is 19.8 ± 9 Gyr. The constraint imposed by stellar popula- 
tions in bright cluster ellipticals (Salpeter IMF) yields an 
age between 12 and 13 Gyr regardless of the prior on ho 
or on the population synthesis model chosen, although im- 
posing a prior on ho results in smaller error bars, roughly 
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around 13.2l^;g Gyr. This is in agreement with estimates 
of the age of the universe from Type la supernovae at high 
redshift (Perlmutter et al. 1999; Riess et al. 1998). Using a 
Scalo IMF does not change the age estimates significantly: 
the best fit gives ages less than 5% compared to a Salpeter 
IMF, i.e. well below the error bars. 

The latest measurement of the age of the oldest stars 
in the Milky Way through the decay of ^''^U gives a value of 
12.5±3 Gyr (Cayrel et al. 2001), which is consistent with the 
above ages if we assume the process of star formation started 
in our galaxy 1 — 3 Gyrs after the Big Bang, corresponding 
to a formation redshift zf ^ 2 for a reasonable range of cos- 
mologies. Another technique which allows for a reasonably 
accurate estimate of the age of our galaxy involves globular 
clusters. The ages of globular clusters can be inferred in a 
distance-independent way by analyzing several features in 
the stellar colour-magnitude diagram such as the luminosity 
gap between the main sequence turnoff and the base of the 
zero age horizontal branch, or the colour gap between the 
turnofi^ point and the tip of the red giant branch. Salaris & 
Weiss (1998) explored a sample of halo and disk globular 
clusters, finding the age of the oldest cluster (NGC 6366) 
to be 12.2 ± 1.1 Gyr. However, Pont et al. (1998) estimate 
an age of 14 Gyr for M92 with a more standard distance- 
dependent technique, involving a fit to the main-sequence in 
a colour-magnitude diagram. M92 is estimated by Salaris & 
Weiss (1998) to be 11.0 ± 1.1 Gyr old. 

The estimates from CMBR data give values that are 
perfectly consistent with the analysis using stellar popula- 
tions. Since the new BOOMERanG data is perfectly con- 
sistent with present estimates of the baryon content from 
Big Bang Nucleosynthesis (Buries et al. 1999) the prior 
flbh'^ = 0.02 ±0.002 has little effect on the results as we can 
see in the fifth row of Table 1. We note that our CMBR anal- 
ysis is restricted to a specific class of models based on adia- 
batic primordial perturbations and with a limited number of 
parameters. Assuming different mechanisms of structure for- 
mation than those predicted by infiation such as topological 
defects and/or isocurvature fluctuations would drastically 
change our conclusions. Furthermore, restricting our analy- 
sis to purely baryonic universes would yield higher values for 
the age of the universe, namely around ^ 22 Gyr (Griffiths, 
Melchiorri & Silk 2001). Thus, the consistency between the 
age values inferred from CMBR and those obtained by stel- 
lar populations can be considered as a further confirmation 
of the standard inflationary scenario. 

As an interesting check of the models, we decided to 
invert the analysis presented here so that the formation red- 
shift and the average metallicity of stellar populations in 
bright ellipticals could be inferred from estimates to the age 
of the universe. With a conservative prior on Bubble's con- 
stant (Ho = 70 ± 15 km Mpc~^), we find the average 
metallicity to be Z/Zq = 1.10 ± 0.18 and 1.17 ± 0.31 for a 
Salpeter and a Scalo IMF, respectively, using the population 
synthesis models of Bruzual & Chariot. Adding the age con- 
straint from CMBR data does not significantly change the 
resuh, giving Z/Zq = 1.14 ±0.18 (Salpeter) and 1.14 ±0.28 
(Scalo). On the other hand, all formation redshifts explored 
in this paper {2 < zf < 10) are allowed, and only the latest 
formation epochs {zf ~ 2) are mildly ruled out at the 68% 
confidence level, which is compatible with the latest esti- 



mates from morphological and spectroscopic studies of high 
redshift clusters (Van Dokkum & Franx 2001). 

Hence, we have presented the study of the colours of 
bright cluster ellipticals as an additional analysis to be in- 
corporated in the medley of cosmological probes. With a 
set of reasonable assumptions for the stellar populations in 
this type of galaxies, we infer an age of the universe of 
13.2^2 Q Gyr and a final result — combining the results 
from stellar populations in cluster ellipticals, the angular 
power spectrum of the CMBR and type la supernovae — 
of 13.4lj^ Q Gyr. Without the assumption of any priors, the 
combined analysis rules out universes older than 18 Gyr at 
a 90% confidence level. 
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TABLE 1: Constraints on the Age of the Universe 

Prior B&C 



Yi 



CMBR CMBR+B&C CMBR+B&C+SN-Ia Database 



h = 0.70 ± 0.15,a;6 = 0.025 ± 0.01 
h = 0.72 ± 0.08,a;b = 0.025 ± 0.01 
h = 0.70 ± 0.15, ivt = 0.025 ± 0.01, U = 
h = 0.70 ± 0.15,aJb = 0.02 ± 0.002 
No prior 
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Figure 1: Age-metallicity diagram for cluster C10024+16 (z = 0.39) using the photometry from Stanford, Eisenhardt & Dickinson 
(1998). The contours are at the 68, 90 and 95%(thick line) confidence levels of the defined in the text. Solid (dotted) lines 
correspond to the population synthesis models of Bruzual & Chariot (2001) and Yi et al. (1999), respectively. Both assume a Salpeter 
initial mass function. The arrows and one of the top axes give the stellar ages for three different formation redshifts assuming a 
A-dominated flat cosmology (ACDM, Qrn = 0.3; ho = 0.7). The other axes on top give formation redshifts for other popular 
cosmologies: OCDM (f2m = 0.3; Qa = 0; /iq = 0.6); and SCDM {Qm = 1; f^A = 0; ho = 0.55), The horizontal line shows the lower limit 

to the metallicity expected for bright cluster ellipticals. 
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Figure 2: Estimate of the age of the universe as a function of Hq preferred by a joint analysis combining the stellar populations in 

bright ellipticals for a Salpeter initial mass function and the observed angular power spectrum of the CMBR. The 238 y 
age-measurement of an old halo star in our galaxy of Cayrel et al. (2001) is also shown as dashed lines, however one should shift this 
age upwards by an amount corresponding to the lapse between 4 = and the first processes of star formation in our galaxy. The age of 
the oldest halo globular cluster in the sample of Salaris & Weiss (1998) is also shown. 
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Figure 3: Likelihood curves for the estimates of the age of the universe using a comparison between the colours of bright cluster 
ellipticals and the latest population synthesis models of Bruzual &; Chariot (B&C; solid line); and the angular power spectrum of the 
CMBR (dashed). The shaded area gives the likelihood of the combined analysis. Universes older than 18 Gyr are ruled out at a 90% 

confidence level. 
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